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1. Introduction 

At the modern level of analytic calculations in elementary particle physics one needs to eval- 
uate thousands and millions of multiloop Feynman integrals. To evaluate a family of Feynman 
integrals which have the same structure of the integrand and differ by powers of propagators (in- 
dices) the standard strategy is to apply integration-by-parts (IBP) relations and solve the problem in 
two steps: a reduction of any given Feynman integral to so-called master integrals and the evalua- 
tion of the master integrals. In this contribution we describe the computation of Feynman integrals 
needed for the evaluation of the three-loop corrections to the static QCD quark potential. In par- 
ticular we describe the use of the Mellin-Barnes technique to evaluate master integrals and the 
application of Grobner bases to solve the IBP relations. 

In Section 2 we present a brief review of the method of Mellin-Barnes (MB) representation 
and exemplify it by the evaluation of a non-trivial integral contributing to the three-loop potential. 
Afterwards we explain in Section 3 the main features of the algorithm called FIRE (Feynman Inte- 
gral REduction) which is based on an extension of the classical Buchberger algorithm to construct 
Grobner bases (see, e.g., Ref. |[l]]). In FIRE Grobner bases are naturally combined with the well- 
known Laporta algorithm. Finally, in Section 4 we present the first results of this evaluation: the 
contributions proportional to «^ and n^, where «/ is the number of massless quarks. 



2. Mellin-Barnes technique 

The MB representation 

{x+yY J-i^ r(A) 271/ ^ ■ ' 

can be applied to replace a sum of two terms raised to some power by their products in some powers. 
For planar diagrams, experience shows that a minimal number of MB integrations is achieved if 
one introduces them loop by loop, i.e. one derives a MB representation for a one-loop subintegral, 
inserts it into a higher two-loop integral, etc. Consider, for example, the dimensionally regularized 
Feynman integral of Fig. |T] which we denote by . . . A straightforward implementation 
of the loop-by-loop strategy leads to a six-fold MB representation which reads 



Figure 1: Feynman integral appearing in the calcula- 
tion of the three-loop corrections to the static potential. 
The straight and wiggled lines correspond to massless 
scalar and static propagators, respectively. The num- 
bers next to the lines refers to the corresponding index 
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(271/) 



. r(a4+zi,2)r(a6+z4,5)r(i/2-z3)r(i/2-z6) 

\X ^ r(a5-Z2,5)r(ai,2,3,4,6.7,8+a9,ll/2 + 2£-4 + Zi,...,6 



r(a9/2 + Z6)r(aii/2 + Z3)r(2 - ^3^4 -an/2-e-zi +Z3)r(aii/2 - 2 + ^3^4.8 + £ + ^1,2,3) 
r(fl;i()/2+ l/2-Z3^6)r(8-ai,...,8 -fl'9,10,ll/2-4£-Zi_4+Z3^6) 
xT{ai^,_^ + fl;9jo,ii/2 - 6 + 3£ + zi,4)r(6 - ai,2.3,4,6,7.8 - ^9,10,11/2 - 3£ - 21,2,4,5) 

X r(2 - a4,8 - ai 1/2 - £ - Z2,3)r(2 - £21,6 - a()/2 - £ - Z4 +Z6)r(2 - flj - aio/2 - £ + Z2,3,5,6) 
Xr(2 - £26,7 - ^9/2 - £ - Z5,6)r(ai,6,7 + ^9/2 - 2 + £ + Z4,5,6) , (2-2) 

where a3,4,8,ii = aj,+aA, + a^ + ci\\, Zi,2,3 = Zi +Z2 + Z3, etc. By definition, any integration con- 
tour over z, should go to the right (left) of poles of Gamma functions with +z-dependence (— z- 
dependence). 

There are two strategies for resolving the singularities in £ in MB integrals suggested in 
Refs. IP (see also Chapter 4 of [Q]). The second one was formulated algorithmically [|5[ |6|], 
and the corresponding public code MB . m has become by now a standard way to evaluate MB 
integrals in an expansion in £. It can be combined with the program AMBRE |^] which can be 
used to derive MB representations in the loop-by-loop approach. Using MB . m and evaluating the 
resulting finite MB integrals by corollaries of Barnes lemmas we have obtained, for example, the 
following result for one of the master integrals occuring in the reduction of F{a\,. . . ,ai\) 



F(l,..., 1,0,1) 



(_^2)3+3e^2 



5671^ 112714 16712,^(3) 8C(5) 

+ ^77^ + + + 0{e) 



135£ 



135 



(2.3) 



m 



At the three- and four-loop level, the method of MB representation was successfully applied 

II, ©nil]- 

The loop-by-loop approach becomes problematic for non-planar diagrams. A typical phe- 
nomenon is that factors such as (—1)"'' arise. This means that the convergence of MB integrals at 
large values of Im(z,) is no longer guaranteed. Moreover, poles in £ arise not only due to "glu- 
ing" of poles of different nature (a typical example is the product r(£ + z)r(— z); there is no space 
between the first poles of the two gamma functions if £ ^ 0) but also from the integration over 
large Im(z,). A safe way to proceed with non-planar diagrams is to start from an alpha (Feynman) 



parametric representation and apply (2.1) to the two basic functions in this representation. (See 



also Ref. [12] for a discussion of this problem.) 



3. FIRE: Feynman Integral REduction 



There are three approaches to solve IBP relations [ ]13| ] in a systematic way: Laporta's algo 



rithmljlj], Baikov's method [ ]15[ ] and two approaches using Grobner bases which are described 
in Refs. [16] and [17, 19, |^], respectively. To solve the reduction problems arising in the 
evaluation of the three-loop potential we use the latter algorithm whose computer implementation 
is called FIRE. It is based on a generalized Buchberger algorithm for constructing Grobner-type 
bases associated with polynomials of shift operators. This method was recently used to evaluate a 



family of nontrivial three-loop Feynman integrals [21]. Let us in the following describe some new 
features of this algorithm. 
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Similarly to Laporta's algorithm, in our approach we work in a given sector, i.e. a domain 
of integer indices a, where some indices are positive and the rest of the indices are non-positive. 
The aim is to express any integral from the sector in terms of master integrals of this sector and 
integrals from lower sectors, where at least one more index is non-positive. It turns out that in the 
higher sectors (with a small number of non-positive indices) the corresponding s-basis [|l8, 2^, 19] 



(a kind of a Grobner basis) can be constructed easily (and, in most cases, even automatically). 

In the opposite situation where a lot of non-positive indices occur, 5-bases are constructed 
not so easily. Usually there is the possibility to explicitly perform an integration over some loop 
momentum for general value of e with results in terms of gamma functions. A straightforward 
way to do this leads to multiple summations and turns out to be impractical. However, there is 
an alternative approach which can be illustrated using the example of diagram of Fig. [l]: consider 
the region a2,a5,aio < and a-],a% > (i.e. the union of the sectors with such restrictions). In 
this situation one can integrate over the middle loop momentum / which enters the propagators 
of the central subgraph with five lines. By constructing an s-basis for this region, it turns out 
possible to solve the IBP relations for the corresponding subintegral over I in order to express any 
such subintegral in terms of master integrals. In other words, the indices a2,a5,aio,aT,ag can be 
reduced to their boundary values, i.e. a2,a5,aio = 0, aT,as, = 1, up to integrals that drop out 
from this region. Then, after using this reduction procedure, it will be sufficient to use explicit 
integration formulae only for the boundary values of the indices. This replacement is very simple, 
without multiple summations. 

It turned out possible to implement the solution of the recursive problem for the subgraph in 
terms of the Feynman integrals for the whole graph. In this reduction, pure powers of the parame- 
ters which are external for the subgraph transform naturally into the corresponding shift operators 
and their inverse. Integrals which are obtained from initial integrals by an expUcit integration over 
a loop momentum in terms of gamma function usually involve a propagator with an analytic regu- 
larization by an amount proportional to e (and, sometimes, 2e). After this integration we obtain a 
two-loop reduction problem with seven indices which is then solved by FIRE. 

After using Grobner bases in higher sectors and an explicit integration in lower sectors, it is 
still necessary to solve the reduction problem in a relatively small number of intermediate sectors. 
In these cases we turn to Laporta's algorithm implemented as part of FIRE. 

4. Evaluating three-loop static quark potential 

The QCD potential between a static quark and its antiquark can be cast in the form 
47iCFas 



(4.1) 



where the renormalization scale of ot^ is set to q^. The one-loop contribution ai is known since 
almost 30 years and also the two-loop term has already been computed end of the nineties. Fur- 
thermore logarithmic contributions are known at three- and four-loop level. Explicit results and the 
references are nicely summarized in the review p^]. The non-logarithmic third-order term, as, is 
still unknown. It is conveniently be parametrized in the form 

as = a^^\] + a'^^n} + a'^^ni + a'^^ , (4.2) 
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where ni denotes the number of massless quarks. Using the techniques described above we evalu- 
ated the coefficients af^ and the CaTp part of a"^^ which read 

J3)_ /20y 3 (2) _/12541 368C(3) 647r4\ 2 .... 
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